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ON THE SMALLER PERTURBATIONS OF THE LUNAR 

ARGUMENTS* 

BY 

ERNEST WILLIAM BEOWN 

The problem of finding the effects of planetary action on the motion of the 
moon has not yet been satisfactorily solved. As far as periodic terms are con- 
cerned the works of Hill, f Newcomb $ and Kadau § have tended to consider- 
ably shorten and simplify the calculation, reducing it, when once the disturbing 
function has been computed, to some simple substitutions in a general formula. 
But the method is not available for finding the non-periodic changes produced 
on the mean motions of the perigee and node. I propose to obtain here some 
simple formulae which permit these perturbations to be found quickly and with 
all the accuracy that can be desired. The results will be applied to the effects 
produced by planetary action and by the figure of the earth, but certain of the 
constants whose values are given below are available for the calculations of any 
non-periodic perturbations of the perigee and node. The notation adopted here 
is, in general, that of my paper " On the formation of the derivatives of the 
lunar coordinates with respect to the elements " || ; the investigation below is, 
in fact, a continuation of that paper to a certain extent. 

(i) Formulae for finding the non-periodic perturbations of the 
perigee and node. 

1. The motion of the moon under the attractions of the earth and the sun 
considered as particles is supposed to have been found. The coordinates are 
thus expressed in trigonometric series of the angles w t , w 5 , w 6 and an angle not 

* Presented to the Society February 27, 1904. Received for publication February 26, 1904. 

f On certain lunar inequalities due to the action of Jupiter, Washington Astronomical 
Papers, vol. 3 (1891), pp. 373-393. 

t Action of the Planets on the Moon, Washington Astronomical Papers, vol. 5 (1894), 
pp. 101-295. 

I Becherches concernant les inigalitSs planetaires du mouvement de la lune, Annales de l'Obser- 
vatoire de Paris (Memoires), vol. 21 (1892), pp. 1-114. 

|| Transactions of the American Mathematical Society, vol. 4 (1903), pp. 234- 
248. This paper will be referred to below by the letter F . 
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containing the constants of integration, with coefficients depending on constants 
of integration c x , c 2 , c 3 and known constants. The angles w. are of the form 

w i = bJ-\-e, w 5 =b s t + Tr, w 6 = b e t + 0, 

where b 4 = n , is the mean motion and b 5 , b e are the mean motions of the perigee 
and node; the b. are functions of the c { and known constants ; e,ir, 6 are the 
other three constants of integration. The c t , w. are called below the " elements " 
of the moon's motion. 

If a new small disturbing function XR be added, the resulting changes in 
lunar elements are obtained by solving the equations (F., art. 18) 

dc. dR dw. dR 

W ^ = X ^ 3 ' ^=- x a^ + ^- <«=i.».»;i=«.«.«)- 

In all practical cases squares of X may be neglected so that the undisturbed 
values of the elements may be substituted in the derivatives of R , — a fact which 
enables us to substitute for R its constant term when R is expanded in a 
trigonometric series, for we are only concerned here with the constant parts of 
dtVj/dt. We therefore use the letter R henceforward to denote the constant 
term of that function. 

The equations (1) can now be solved and they give 

dR T 

c. = const., w.= — Xt~^ 1- b.t + const. 

' J dc. ' ' 

1— 3 

The theory is at first supposed to be expressed, not in terms of w., c ; , but in 
terms of w and n, c 2 , c 3 , the semi-canonical system (F., art. 10): c 2 contains 
the square of the eccentricity constant as a factor, and c 3 that of the inclination 
constant. Let Bn, Bc 2 , 8c 3 be the changes in n, c 2 , c 3 caused by Xi? and Bb f 
the changes in the b . Then, since 6 4 = n, 

8b. = — X^^- + hn. 

The old definition of the letter n is to be retained, that is, it is to represent 
the mean motion ; thus Bb t = and Bn = XdR/d^ . We shall, for the present, 
define c 2 , c 3 to be the same as when X = , so that Bc 2 = , Bc 3 = . These 



1904] PERTURBATIONS OF THE LUNAR ARGUMENTS 281 

definitions are possible because c 15 c 2 , c 3 are constant whether X is zero or not. 
Hence 



».. 



_ fdR db 5 DR\ __ /cLS dV dR\ 

• ~ ~ \ dc 2 ~ dn 3c, / ' 6 — — \ 3c 3 ~ dn dc e J ' 

where Z> 5 , 5 6 are expressed in terms of n, c 2 , c 3 and i? is expressed in terms of 



I have shown (F., no. 10) that 

dn dc 2 ' dn dc B ' 

where b 2 , c, are expressed in terms of n, c 2 , c 3 . Hence 

(2) ^="^lf' S& a=-^|f. 

where ^ is now expressed in terms of n c 2 , c s . 
These are the required expressions for 8Z> 5 , SZ> 6 . 

3. The results (2) are remarkable in that they show the greater simplicity 
obtained by using the semi-canonical system n, c 2 , c 3 instead of the canonical 
system c, , c 2 , e 3 ; a similar result was found in another connection (F., art. 10). 
Moreover, the main difficulty in obtaining accuracy, namely, the formation of 
derivatives with respect to n, has disappeared. One or two transformations 
are all that are necessary in order to calculate Sb 5 , Bb 6 with good accuracy from 
a numerical theory and with complete accuracy from a theory in which the 
numerical value of the ratio of the mean motions is substituted while the other 
constants are left in a literal form. The results of my theory are therefore 
immediately available when the derivatives of It with respect to c 2 , c 3 have 
been expressed in terms of derivatives with respect to e , k — the constants of the 
theory. 

In the applications it is never necessary to go beyond the second powers of 
the eccentricities and inclinations, and all the functions considered are expan- 
sible in squares of these quantities and of the ratio of the parallaxes ; the 
square of this latter ratio is negligible. As we have to form derivatives with 
respect to c 2 , c 3 , that is, with respect to e\ k 2 , it will be necessary, for this 
degree of accuracy, to find the expansions as far as the second powers and 
products of e 2 , e' , k 2 . 

4. We have, when the independent variables 6,, c 2 , c 3 are changed to 5 X , 
e 2 k 2 

a? ~ ~dc~ 2 &T 2 + e^ ae 2 ' ak 2 ~ dc 2 ek 2 + ac 3 3k 2 ' 

where d/de 2 , d/dk 2 denote derivatives with respect to e 2 , k 2 . Since c 2 contains 
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the factor e 2 , and c 3 the factor k 2 , we obtain from (2), on omitting the terms 
of an order higher than the second, 

£5 Si? 3c, 
Sb. dR Be 2 3k 2 3e 2 

5 



X 3c 2 dc 2 8c 2 dc 3 
3e 2 Be 2 3k 2 
with a similar equation for Bb 6 . 

When R, c 2 , c 3 are expanded in powers of e 2 , e' , k 2 , in the form 

R = n 2 a 2 [-S(O) + e 2 R(e 2 ) + k 2 i?(k 2 ) + •••]» 
c 2 = na 2 [e 2 c 2 (e 2 ) + e 4 c 2 (e 4 )+ eV 2 c 2 (e 2 e' 2 ) + •••], 
c 3 = na 2 [k*c 3 (k 2 ) + k 4 c 3 (k<)+ e 2 k 2 c 3 (e 2 k 2 ) +•••], 

where i?(0), • • •, c 2 (e 2 ), • • •, c 3 (k 2 ), • • • are functions of n' /n only and there- 
fore numerical coefficients, we obtain finally 

S6 5 -S(e 2 )(" 2e 2 c 2 (e 4 ) + k 2 c 2 (e 2 k 2 ) + e' 2 c 2 (e 2 e' 2 )-j 



b- 



\u c 2 (e 2 ) L c 2 (e 2 ) J 

^ ig(k 2 )c 3 (e 2 k 2 ) 2e 2 -R(e 4 ) + k 2 i?(e 2 k 2 ) + e'lg(eV 2 ) 

" c 2 (e 2 )c 3 (k 2 ) + ~ c 2 (e 2 ) ' 

with a similar expression for Bb 6 in which c 3 , c 2 , k, e are respectively put for 

These expressions are not completely accurate to the orders given owing to 
the final definitions of e, k. These constants are defined by means of the prin- 
cipal elliptic inequality and the principal term in latitude respectively. In gen- 
eral the periodic terms in R will produce in the coefficients of these terms addi- 
tional parts of the form Xa e, X/3 k where a , /3 are numerical coefficients. In 
order to reduce to the final definitions it is therefore necessary to replace e by 
(1 — Xa )e, k by (1 — X/3 )k. The change in Sb 2 is therefore 

-2^§X« e 2 -2||x/3 k 2 . 

As Bb b jBe 2 , 36 5 /3k 2 are of the order (n' fn) 2 at least, the orders of these addi- 
tional parts relatively to the corresponding parts of order e 2 , k 2 already found, 
are in the ratio (n' /n) 2 : 1 . A similar remark holds in the case of Sb e . In no 
actual case do these portions give additions to the annual mean motions so great 
as 0".005 ; they therefore can and will be neglected. These additional parts, 
however, do not violate the general theorem of this paper, namely, that algebra- 
ical developments in powers of n' jn are not required, for the coefficients a , /3 a 
can be obtained without such developments (F) . 
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5. The following numerical values are needed : 

\ = - .11845 e 2 - .02324 e 4 - .26363 e 2 k 2 - .00110 e 2 e'\ 
nar 

\=- 2.00206 k 2 - 1.96376 k 4 - .28546 e 2 k 2 - .00568 e' 2 k 2 ; 
na l 

n = 17325594"; 

e = . 10983, k=. 044892, e' = .01677. 

The values of c 2 , c 3 have been found by the method given in F., art. 11. 

The relations between my constants and those (e, 7) of Delaunay's theory, 
though not needed here, may be given : 

e = + 2.00054 e - .36682 e 8 - 2.01160 e-y 2 + .04873 ee'\ 

k = + 1.00013 7 - .49609 7 s - .49924 ye 2 - .00037 ye' 2 . 

(ii) Application to planetary action. 

6. Let x, y, z be the coordinates of the moon, £, 77, f those of a planet re- 
ferred to moving axes of a:, y in the plane of the ecliptic of which the aj-axis 
points to the mean place of the sun, and an axis perpendicular to this plane. 
Then if r 2 = x 2 -f y 2 -f z 2 , D 2 = I 2 -f if + f 2 and \ denotes the mass of the 
planet, we have, neglecting the ratio of the parallax of the sun to that of the moon, 

„ _ 3 («£ + yv + zSf 1 f_ 

u — 2 jy> 2 jys' 

It will tend to simplicity if we neglect at the outset certain quantities whose 
squares are certainly negligible in the results ; these are the solar eccentricity 
and the product of the inclinations of the planet and the moon to the ecliptic. 
As the constant term in S, is alone required we can further neglect the products 
xy^ij/IP, xztjZ/IP, yzrjZ/D 5 . The function then can be put into the form : 

£ = i { x 2 + y 2 -2z 2 )(~-^ + Ux 2 -y 2 ) 1 ^ 

= H^-2* 2 ) (^ 3 -^) o + tK + * 2 )o(^ 2 ) o 

where u = x + y V — 1 , s = x — y V — 1 , and ( ) denotes the constant term 
in the expansion of the included function. The separation of the parts depend- 
ing on the coordinates of the moon from those depending on the coordinates of 
the planet is possible because when the solar eccentricity is neglected, they have 
no common angles. 
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Let n', a', be the mean motion and mean distanoe of the sun ; w, a, those 
of the moon ; m', the mass of the sun ; m = »'/( n — n')\ a, a certain constant 
which differs little from a; 

then, since (« 2 ) = (s 2 ) , the expression for J? may be written 

m 2 a 2 , ,r /ws-2a 2 \ , /« ! \1 

In this expression the only parts containing e, k are the coefficients of H, K 
within the square brackets. I already have the results for these since they were 
required in obtaining the motion in the main problem. I find 

( ^~a" 2 ~ ) + - 38229 e * ~ 6 - 00851 k ' + - 0398 e * 

- 3.0395 e 2 k 2 + .0053 k 4 + ■••, 

• + .25496 e 2 + .18030 k 2 - .0142 e 4 

-. 4094 e 2 k 2 -. 2035 k 4 +-.-. 



©.- 



7. We must next calculate IT, K. First, owing to the numerical values of 
the eccentricities, inclinations and masses of the planets, the only possible sensible 
effect of the eccentricities and inclinations is the inclination of Venus. Let 
then V", V be the mean longitudes of a planet and the sun, a" , a their mean 
distances, y" the sine of half inclination of the planet's orbit to the elliptic, h" 
the longitude of the planet's node. Then, omitting powers of 7" beyond the 
second, 

I = - a' + (1 - 7 " 2 )a" cos ( V" - V) + j" 2 a" cos ( V" + V - 2h"), 

v = (l-j" 2 ) a" sin (V"-V)- y" 2 a" sin V" + V - 2h" ) , 

£=2y"a" sin(V"-h"); 

B 1 = a' 2 + a" 2 - 2a a" (1 - y" 2 ) cos( V" - V) - 2y" 2 a'a"cos ( V" + V'-2h") 

= B\- 2y" 2 a'a"( V" + V- 2h"), 

?- v 2 = a' 2 -2a'a"(l - y" 2 ) cos (V- V) + a"[{l- 2 7 " 2 )cos2( V" - V) 

- 2y" 2 a" { a' cos ( V" + V- 2h") - a" cos 2{V- h")} f 

^ = 2 7 " 2 a" 2 { 1 - cos 2 ( V" - h" ) } . 



1904] 
Let 

Then 



PERTURBATIONS OF THE LUNAR ARGUMENTS 

D- n a' H = i« + ZAf cos ( V" - V) + 2Af cos 2 ( F" - V ) + 
The following table gives the values of the various quantities : 
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Planet. 


Mercury. 


Venus. 


Mars. 


Jupiter. 


Saturn. 


T 


6000000 


408000 


3093500 


1047.35 


3501.6 


, a" 
log ~y 


9.5878 


9.85934 


.1829 


.71623 


.9795 


A™ 


1.4359 


4.9317 


1.0858 


10-» X 7.7267 


10- 3 X 1-181 


A? 


2.5658 


41.9340 


2.3140 


10-* X 3.3035 


10- 5 X 1-356 


Af 


1.9560 


40.7575 


2.2042 


10-* X 1.4856 


10- 6 X 3-483 


Af 


1.2086 


37.8764 


1.9597 


10- 5 X 4.8834 


10-' X 6.35 



Of these numbers, the values of the masses are taken from the best modern 
determinations ; the ratios of the mean distances are certain to the degree of 
accuracy required here ; the values of A ( ^ for Venus and of Af } for Mars have 
been calculated from the tables given by Newcomb ; * the rest of the numbers 
are taken from the corresponding table of Hansen, f 

It is directly seen that the inclination of Venus exerts a slight effect. For 
since ?" = .0296, the ratio of the second term of H to the first is 1/43, and the 
part produced by the first term is about 1" as we shall see in the following table. 
The effect on K is smaller. 



* Action of the Planets on the Moon, Washington Astronomical Papers, vol. 5 (1894), 
pp. 248-261. 

f Darlegung der theoretischen Berechnung der in den Mondtafeln angewandten Stbrungen ( Abh. 1), 
Abliandlungen der mathematisch-physischen Classe der Kdniglioh Saoh- 
sisohenGesellschaftederWissenschaften, vol. 6 (1862), p. 487. The table above is not 
quite consistent for the values Af 1 for Venus and that of A^~> for Mars include the squares of the 
inclinations and eccentricities of those planets, while the other ^<™> do not. But as all these 
squares are negligible, except that of the inclination of Venus, the slight inconsistency makes 
no difference in the results. 

Trans. Am. Math. Soc. 19 
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8. The results of the calculations are as follows : 
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Planet. 


W>H 


Part from H. 


1W 


Pari from K. 


865 


S&6 


«&6 


«>6 


Mercury, 

Venus, 

Mars, 

Jupiter, 

Saturn, 


+ 24 

1181 

35 

738 

34 


+ 0".018 
.906 
.027 
.566 
.026 


— 0".017 
.860 
.025 
.537 
.025 


+ 20.6 

687.9 

4.7 

10.2 

0.0 


+ 0".031 

1.060 

.007 

.016 

.000 


+ 0".001 
.041 
.000 
.001 
.000 


Sums, 


4-2012 


+ l // .543 


— 1".464 


+ 723.4 


+ 1 // .114 


+ 0".043 



The complete annual mean motions contributed by direct planetary action 

are therefore 

S6 6 = + 2".66 , Bb 6 1".42 . 

In order to show the degree of accuracy obtained by the method, I give the 
parts of these contributed by e 2 , k 2 . They are : 

in S6 5 , from H, by e 2 , - 0".004 ; by k 2 , - 0".025 ; 

» " " K, " — 0".005; " — 0".009; 

"S& 6 , " H, " -0".009; » +0".006; 

K, " -0".003; » .000. 



u 11 u 



The parts independent of e 2 , k 2 have also been calculated directly by the method 
of my lunar theory, and the results agree exactly, as they should, with those 
obtained above. 

(iii) Application to the action arising from the figure of the earth. 

9. The disturbing function has been fully expanded by Hill.* The only 
terms which will produce constant portions are, in his notation, 

§H!-2 7 2 + 2 7 *]J 3 +2^[ 7 2 - 7 *]pcos2», 



\ = 



£ 



or, in rectangular coordinates, 

where r 2 = x 2 + y 2 + z 2 . In this expression X is a small coefficient determined 
from observation. I shall here adopt Hill's value f for it, since for any other 

* G. W. Hill, Determination of the inequalities in the moon's motion which are produced by the 
figure of the earth, Washington Astronomical Papers, vol. 3 (1891), p. 213. 
t Ibid., p. 340. 
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value the changes to be made are quite simple ; the most probable value for X 
is discussed elsewhere. * 

From the results of my lunar theory, I find 

/ 3 2 3 \ 

( |_ _ *.JL \ = . . . + .1287 e 2 - 2.0060 k 2 + .0571 e* - 1.0817 e 2 k 2 

+ .0448 k 4 + .0407 eV 2 - .0136 kV 2 + • • •. 
With X = 10~ 7 x 3.6864, the results, on substituting in (3), are 



Characteristic. 


1 


e2 


k^ 


e« 


Sums. 


<% 5 
<56 6 


+ 6*. 936 
— 6*400 


+ 0".041 
— O".042 


— 0".118 
+ O".02.'i 


+ 0".001 
.000 


+ e'.seo 
— &>a\g 



The values for S6 5 , Bb s differ but little from those of Hill,| + 6".820, and 
— 6".413, obtained by Delaunay's method with a literal expansion in powers 
n' fn. The terms of characteristic unity have also been obtained by using the 
method of my theory, and the results agree with those just found. $ 

Havekford College, 
February 23, 1904. 



* E. W. BROWN, On the degree of accuracy of the new lunar theory and on the final values of the 
mean motions of the perigee and node, Monthly Notices of the Royal Astronomical 
Society, vol. 64 (1904), pp. 524-534. 

t Ibid., p. 344. 

% The complete results for the theoretical motions of the perigee and node, together with a 
comparison with observation, will be found in the paper On the degree, etc., to which reference 
has just been made. 



